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INTRODUCTION 


For  the  fundamental  two-receiver  passive  location  system  illustrated 
in  figure  1,  the  available  time-delayed  signals  are  observed  in  the  presence 
of  receiver  noise  and/or  interference.  By  temporarily  neglecting  possible 
doppler  effects,  the  received  signals  can  be  expressed  in  the  simple, 
general  form 


r  L  (t)  =  s(t)  +  ^(t) 

0  <:  t  <  T 

(1) 

r,2(t)  =  As(t  -  r  )  +  n2(t) 

T  s  t  <  T.  +  T 
o  2  o 

(2) 

where  s(t)  is  an  unknown  signal,  \  is  an  unknown  relative  amplitude, 
r  =  (R„-  R..)/c  is  an  unknown  relative  time-delay,  and  the  n  (t)  are  in- 

u  \  1 

dependent  noise  and/or  interference  present  at  each  physically  separate, 
time- synchronous  receiver.  The  complicating  effects  of  signal  doppler, 
which  may  or  may  not  be  present,  are  reviewed  in  a  later  section. 

This  paper  considers  the  case  of  bandpass  signals  with  one-sided 
spectral  bandwidths  of  W  Hertz  centered  about  the  angular  carrier  frequency 
u  radians/sec.  (Similar  results  can  be  obtained  for  lowpass  signals. ) 
Consequently,  the  observed  signals  can  alternatively  be  expressed  in  terms 
of  their  complex  (phasor)  amplitude  as 

?1(t)  =  s(t)  +  n1(t)  =  I1(t)+  JQ^t)  0  <  t  s  Tl  (3) 
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Figure  1.  Time-Delayed  Signals  Observed  in  the  Presence  of  Noise  at  Two  Locations 


?2(t)  =  xi(t-r)ej^  +  n2(t)  =  I2(t)  +  JQ2(t)  tq  s  t  x  Tg  +  tq  (4) 

where  the  overbars  denote  complex  variables;  e^  =e  ^WoT  is  an  unknown, 
ambiguous  relative  phase  distributed  uniformly  over  [  0,  2 7r  ]  radians;  and 
the  I  (t)  and  Q^(t)  are  the  complex  quadrature  components  observed  at 
oach  receiver. 

Since  digital  methods  of  transferring  and  processing  the  received 
signals  are  of  interest,  discrete-time  samples  of  the  form 

1=1,2 . Nx 

(5) 

P  =  1,  2, .  .  .  ,  Ng 

(6) 

where  T  =  N  At  and  T  =  N  At  are  assumed.  As  long  as  At  <  l/W, 

XX  &  Ld 

the  sampling  theorem  for  time-limited,  bandpass  signals  [  1]  indicates  that 
all  the  information  observed  at  the  receivers  can  be  specified  by  the  com¬ 
plex  sample  sequences 


■u  -  ri(Lit) '  8i +  h.i  +  )Qi,i 


r.  =r0(pAt  +  T  )  =  xi  ^  e^+  n  =1  +  jQ 

2,p  2'  o  /  p+m  2,p  2,p  2, 


■  IN  - 
'^ii-r*1 


d  j-  |N2 

and  r  =  r_ 

I  2,pi  p  =  i  2 
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where  the  underbars  denote  vector  variables.  As  a  simplifying  assumption, 

It  is  assumed  in  equation  6  that  the  unknown  offset  lead  (negative  lag)  Index 

m  =  (t  -  r  )/At  =  -  (t  -  tq ) /At  (7) 

is  an  Integer.  This  assumption  keeps  the  analysis  simple  and  general  and 
can  be  properly  accounted  for  by  a  statistical  processing  loss  after  a  specific 
receiver  bandpass  characteristic  and  a  specific  sampling  Interval,  At,  have 
been  chosen.  For  Nyquist  sampling  rates.  At  =  l/W,  this  loss  will  be  on 
the  order  of  1  dB  [  2  ] .  If  desired  the  loss  can  be  reduced  as  low  as  desired 
by  oversampling  or  can  be  eliminated  by  interpolative  processing  since  the 
two  complex  sample  sequences  specify  all  the  observed  information. 


OPTIMAL  WHITE,  GAUSSIAN  SYSTEM 


If  the  signal  and  noises  are  stationary,  white,  zero-mean,  Gaussian 
stochastic  processes,  the  optimal  likelihood  ratio  statistic  to  detect  the 
presence  of  a  correlated  signal  (hypothesis  2;  -T  £  r-  tq  <  T  )  versus  the 
absence  of  a  signal  (hypothesis  0;  s(t)  =  0)  or  the  presence  of  an  uncorre¬ 
lated  signal  (hypothesis  lj  r-  tq  <  -T  or  t-  Tq  >  T  )  is  [  3]  the  magni¬ 
tude  (envelope)  of  the  normalized  sample  cross-correlation  function  between 

the  two  separately  observed  complex  sequences,  )R-  -  (m)|.  The 

rlr2 

normalized  sample  cross-correlation  function  for  a  candidate  lag,  k,  is 


R=  ;  <k)  = 

fir2 


N(k)  +  k 

N(k)  2  ri,ir2,i-l 


0  £  k  s  (Nx-  1 


i  =  k+  l 


N(k) 

N(kj  S  ri,ir2,i-l 


■K-1) 


k  £  0 


where 


N(k)  = 


min  1 

(N2'Nf  k) 

VI 

VI 

o 

min 

(n  N  -  |k|) 

tH 

1 

1 

\  1,  2  1  '/ 

\  2  / 
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is  the  number  of  overlapping  complex  samples  in  the  sequences  for  a  lag  of 

k  and  the  asterisk  denotes  complex  conjugation.  Note  that  the  maximum 

value  of  N(k)  is  min  (N  ,N  )  and  is  maintained  over  an  interval  of  length 

(IN  -  N  |  +  i)At.  The  optimal  two-receiver  statistic, 

1  ^ 

R-  -  (k)  b  R(k)  , 

1  2 


can  test  for  the  presence  of  signals  with  time-delays  in  the  interval 
[  -T  st-  t  s  T  ]  and  is  strictly  optimal  only  when  N(k),  the  number 

1  O  2 

of  independent  complex  terms  in  the  cross-correlation  sum,  is  large  in  the 
law  of  large  numbers  sense.  Not  only  does  the  cross-correlation  statistic 
lose  its  condition  of  optimality'  near  the  interval's  edges,  but  its  performance 
degrades  unacceptably  since  the  cross-correlator  output  signal-to-noise 
ratio,  r,  for  each  candidate  lag  of  k  is  [  3  ] 


where  T(k)  =  N(k)  At  is  the  overlap  time  for  the  lag  of  k  and 


E  [|i(t)|  2 
E  [|n1(t)|2] 


(ID 
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and 


E  ^ {  Xb  (t  -  r)}2  E  T  |  x  s  (t  -  r)|  2 

^  E[”2<t)]  -[M*] 

are  the  Input  signal-to-noise  ratios  at  receiver  1  and  receiver  2.  The 
probability  of  detection  and  the  probability  of  false  alarm  are  interrelated 
as  [3] 


Pd(k)  =  Q  j  y^nk)  ,  J- 2  in  [pf&(k)j  j  (13) 

where  Q(o,0)  is  Marcum's  Q  Function  {4]. 

When  the  presence  of  a  time-delayed  signal  is  indicated,  the  maximum 
likelihood  estimate  of  the  time-delay  corresponds  to  the  lag,  k,  at  which 
the  interpolated  magnitude  of  the  normalized  sample  cross-correlation  func¬ 
tion  is  maximum  and  is  obtained  from  [  3  ] 


t  -  t  -  mAt  .  (14) 

o 

When  the  output  signal-to-nolse  ratio  is  large,  the  time-delay  estimation 
error  is  asymptotically  unbiased,  Gaussian,  and  efficient  with  a  variance 
of  [3] 


oN - (15) 

T  T  (m)0 

where  0  is  the  rras  bandwidth  of  the  receivers  (about  the  carrier)  in 
radians/second. 
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Figures  2  and  3  illustrate  the  optimal  two-receiver  coherent  passive 
location  system  as  developed  above.  Because  of  the  unknown  signal  and 
noise/interference  levels,  an  adaptive  threshold,  Tq,  must  be  generated 
for  each  candidate  lag  k  as  illustrated  in  figure  3.  Note  that,  as  shown  by 
equations  10,  13,  and  15,  the  detection  and  tim^-delay  estimation  perform¬ 
ance  of  the  optimal  system  at  each  possible  lag,  k,  is  completely  character¬ 
ized  by  the  output  signal-to-noise  ratio,  T(k),  which  increases  linearly  with 
the  overlap  time-bandwidth  product.  Equation  10  indicates  that  the  output 
signal-to-n^ise  ratio  observed  at  the  cross-correlator  will  always  be  less 
than  the  overlap  time-bandwidth  product. 
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Figure  3.  Optimal  Detection  and  Time- Delay  Estimation 


GENERAL  OPTIMALITY 


If  the  signal  and/or  noise  statistics  are  non- Gaussian  but  stationary 
and  N(k)  Is  large,  the  normalised  sample  cross-correlation  function,  R(k), 
will  remain  optimal  in  a  least  mean  square  sense  and  will  attain  the  perform¬ 
ance  of  the  optimal  Gaussian  case.  If,  in  addition,  the  signal  and/or  noise 
statistics  are  not  stationary  (e.  g. ,  pulsed  or  otherwise  deterministically 
amplitude  or  phase  modulated),  the  normalized  sample  cross-correlation 
function  will  still  remain  optimal  in  a  least  mean  square  sense  and  will  still 
attain  the  performance  of  the  optimal  Gaussian,  stationary  case  if  the 
overlap  time,  T(k)  =  N(k)At,  is  sufficiently  long  so  the  mean  square  sample 
statistics  of  the  signal  and  noises. 


and 


N(k) 

N(kj  23 

i=l 


2 


N(k) 

1  1 2 
N(k) 

1=1 


1 

N(k) 


f 
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are  essentially  stationary  (independent  of  the  Interval  length;  i.  e. ,  constant) 
and  the  sample  cross-correlations  between  the  signal  and  noises  are 
essentially  aero. 

The  previous  discussion  indicates  that  a  coherent  passive  location 
system  that  cross-oorrelates  signal  segments  observed  at  physically  sepa¬ 
rate  receivers  is  inherently  robust  (e.  g. ,  will  perform  well  against  either 
pulsed  or  continuous  signals  which  may  be  either  random  or  deterministic) 
and  makes  optimal  use  of  the  available  Information  as  long  as  the  observa¬ 
tion  intervals  and  cross-correlation  intervals  are  statistically  long.  To 
obtain  this  optimal  performance  and  versatility,  the  two- receiver  system 
can  be  implemented  with  either  collocated  or  remote  processing  as  shown 
in  figures  4  and  5. 
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DOPPLER  COMPENSATION 

The  previous  analysis  assumes  that  doppler  effects,  due  to  relative 
motions  between  the  signal  source  and  either  or  both  receivers,  are 
negligible.  Doppler  effects  can  be  accounted  for  by  modification  of 
equations  4  and  6  to 

to  jwd(t_T)  - 

r0<t)  =  \s(t-  r)  e  +  n  (t)  (] 


to  _ 

r„  =  As  e v  e  + 

2,p  p  +  m  2,  p 


where  the  relative  (differential)  doppler  frequency  between  the  two  time- 
delayed  signals  in  radians/second  is 


-“o  (K  -  K ) 


■  v  -  v 


This  simple  formulation  assumes  that  the  signal  bandwidth  W  is  large 

compared  to  each  received  doppler  frequency,  u>.  and  cu  ,  and  that 

dl  d2 

signal  compression/expansion  effects  and  effects  due  to  higher-order 
motion  terms  (e.  g. ,  acceleration)  are  negligible  over  the  prospective 
cross-correlation  overlap  times. 

When  no  signal  is  present  (hypothesis  0)  or  when  the  observed  signals 
are  uncorrelated  (hypothesis  1),  the  linearly  increasing  doppler  phase 
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introduced  into  the  normalized  sample  cross-oorrelation  function  statistic, 
R(k)  (see  equations  5,  8,  and  17),  will  not  cause  a  change  in  the  cross-cor¬ 
relator  output  signal- to -noise  ratio  (equation  10)  since  the  phases  are  already 
random.  For  the  correlated  signal  case  (hypothesis  2)  the  loss  in  cross- 
correlator  output  signal-to-noise  ratio  for  a  particular  doppler  frequency,  u)j, 
and  a  prospective  time-delay  lag,  k,  can,  by  a  modification  of  equation  10, 
be  expressed  as 


E 

|R(k)|2/H2,  wd 

- 

E 

R(k)|  2/Hx  ] 

E  [ 

|R(k)|2/H2,  wd  = 

J 

-  E 

[|R(k)|2/H1] 

(19) 


By  an  examination  of  the  components  of  equation  19,  the  doppler  decorrelation 
loss  factor,  0  s  <  1,  where  -  1  is  the  case  of  no  loss,  can  be  shown 
to  satisfy  the  bound 

n 

_  2 

>j!  e 


N(k) 


U  1  = 

1 

1  J 

N(k) 

2* 

E 

E  '-v2 

- 

i=  1 

N(k)  N(k) 


l;i|2|5c 


-ju;dAt(i-q) 

e 


i=l  qel 


N(k)  N(k) 

_ 

•"t1 

2  ,  -  ,2 

flq 

i  =  l  q  = 1 

.  21 

2 

sin  ^N(k)  YAt) 

N(k)  6 

lUA  \  ! 

i=l 

N(k)  sin  At  j 

(20) 


i 
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where  the  last  Inequality  holds  for  moderate  or  large  N(k)  and  is  proved 
in  the  appendix.  This  doppler  decorrelation  loss  bound  is  illustrated  in 
figure  6.  Note  that  3  dB  of  loss  is  obtained  for 


u . 


If  -  — 

1  d1  2tt 


0.  443 
N(k)  At 


(21) 


If  expected  doppler  decorrelation  losses  are  unacceptable,  they  can 
be  reduced  as  low  as  desired  by  a  doppler  compensation  filter  bank.  Each 
doppler  filter  bank  section  compensates  for  a  prospective  doppler,  a>c,  and 
implements  a  compensated  normalized  sample  cross-correlation  function 
of  the  form 


R(k)  = 


N(k)+k 

i  y  -  -* 

N(k)  Z-i  ri,ir2,i-l 
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which  is  a  Fourier  Transform  of  the  N(k)  normalized  cross-correlation 

_  _* 

product  samples  for  a  lag  of  k,  r^  j  and  can  be  implemented  for 

a  discrete  number  of  equally  spaced  compensated  frequencies 


u> 


A  /  2* 

r4“  =  r  NS 


r=  1,  2,  N 


where 


N  =  min 
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(23) 
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Doppler  Decorrelation  Loss  Bound 


by  two  complex  shift  registers  of  length  N,  N  complex  multipliers,  and 
an  N  Point  Discrete  Fourier  Transform  as  Illustrated  in  figure  7. 

Doppler  compensation  processing  will  reduce  output  signal-to-noise 
ratio  loss  to  L^(k,  )  where  w  is  the  closest  compensation 

frequency  to  If  the  differential  doppler  is  uniformly  distributed  over 
the  compensated  frequency  band,  the  average  doppler  decorrelation  Iobs 
factor  for  a  lag  of  k  will  be  bounded  by 
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gure  7.  Doppler  Compensation  Processing 


DATA  TRANSFER  REDUCTION  BOUNDS 


Direct  implementation  of  the  optimal  two- receiver  coherent  passive 
location  system  developed  above  vill  require  the  transfer  of  one  or  both  of 
the  complex  sequences,  r  and  iv ,  with  a  high  degree  of  accuracy  (large 
number  of  bits  representing  each  receiver  sample).  Since  higher  place  bits 
carry  diminishing  quantities  of  information,  it  is  clear  that  truncation  of 
high  place  bits  will  reduce  data  transfer  requirements  without  significantly 
reducing  system  performance.  Techniques  for  minimizing  the  data  transfer 
required  to  maintain  reliable  system  performance  are  desired.  The  general 
one-channel  data  transfer  reduction  situation  is  illustrated  in  figure  8. 

Rate  distortion  theory  provides  bounds  on  data  transfer  reduction  for 

particular  types  of  Information  sources  and  distortion  (error)  measures. 

The  information  of  a  memory  less,  Gaussian-di  stributed  source  with  a  power 
2 

of  a  car.  be  transferred  over  a  noiseless  channel,  as  illustrated  in  figure  9, 
with  a  mean  square  error  between  the  input  and  output  complex  sequences  of 


D  =  E 


[>  -  .2' 
1  1  i 


(26) 


by  the  use  of  an  information  rate,  R,  of  [  5  ) 


R  =  -log. 


0  £  D  s 


(27) 


\o  / 

where  the  rate  is  in  bits  per  complex  sample.  This  rate  is  also  sufficient 
for  any  non- Gaussian  source  (  5] . 
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One-Channel 


Figure  9.  Rate  Distortion  Theory  Model 
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For  the  one-channel  data  transfer  reduction  situation  of  figure  8  where 
the  data  of  receiver  2  Is  assumed  to  be  quantized,  the  cross-correlator 
output  algnal-to-noise  ratio  (see  equations  8  and  10)  can  be  expressed  In 
the  form 


where  q  is  the  complex  quantizing  noise  of  the  r  sample.  The 
2»p  2,p 

loss  In  cross-correlator  output  signal-to-noise  ratio,  0  ^  s  1,  for  a 
particular  lag,  k,  can  be  shown  to  be  of  the  form 


L  (k)  = - ~ 

a  2 
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and  D/cr  is  the  ratio  of  mean  quantising  noise  power  to  mean  quantiser 
input  signal  power.  The  inequality  of  equation  30  holds  due  to  the  known 
nonpositive  correlation  between  quantiser  input  signal  and  quantising  noise 
when  the  quantizer  bin  outputs  are  chosen  as  the  mean  input  values  of  each 
bin  [6]  or  when  other  more  general  conditions  are  satisfied.  Thus,  the 
quantization  loss  factor  has  been  shown  to  be  bounded  by 


L 

q 


1  +  2 


-R 


(32) 


where  the  last  inequality  is  obtained  by  the  application  of  the  bound  of 
equation  27.  The  rate  distortion  bound  of  equation  32,  which  holds  for  all 
input  signal-to-noise  ratios,  is  shown  in  figure  10. 

If  the  -eceivod  signals  and  noises  are  Gaussian,  performance  within 
1/2  bit  per  complex  sample  of  the  rate  distortion  bound,  as  shown  in 
figure  10,  can  be  obtained  for  all  input  signal-to-noise  ratios  by  either  [  5  ] : 

1.  Optimal  but  complex  entropy  encoding  combined  with  the 
use  of  uniform  I  andQ  quantization  bins,  or 

2.  Uncoded  (direct)  transfer  combined  with  the  use  of  optimal 
minimum  mean  square  Max  (  7 )  quantization  bins. 

Either  of  these  techniques  would  be  moderately  complex  to  implement  and 
both  assume  knowledge  of  the  received  signal  and  noise  probability  densities. 
Performance  would  degrade  significantly  unless  adaptive  processing  were 
also  provided  to  adjust  for  unknown  amplitude  and/or  distribution. 
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CANDIDATE  DATA  TRANSFER  REDUCTION  TECHNIQUES 


The  performance  of  two  simple,  robust  data  transfer  reduction 
techniques  was  obtained  by  digital  simulation  for  the  one-channel  situation 
of  figure  8  with  stationary,  white,  Gaussian  signal  and  noise  sequences. 

Both  techniques  utilized  direct  decimal-to-binary  encoding  (uncoded  transfer). 


The  first  data  transfer  reduction  technique  transferred  and  cross- 
correlated  only  the  quantised  phase  of  each  complex  sample  of  receiver  2. 
The  simulated  phase  cross-correlation  was  of  the  form 
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where  the  tilde  denotes  a  quantized  version  of  the  corresponding  received 
phase  samples 


,  -1 
tan 

and  9n  =  tan  * 
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Uniform  phase  quantization  bins,  as  illustrated  in  figure  11,  were  utilized 

with  the  center  of  each  bln  used  as  the  quantizer  output.  Because  of  the 

assumed  uniformly  distributed  random  receiver  phase,  0^  ,  direct  encoding 

with  uniform  bins  achieves  the  optimal  rate  distortion  bound  for  Gaussian 

2 

Thus,  phase  can  be 


sources  of  R  =  -  -  log^ 


.  2  rr 

where  cx0  =  — 


transferred  with  minimal  distortion  by  extremely  simple  quantization  and 
encoding  techniques.  Simulation  results  obtained  for  the  case  of  0  dB  input 
signal-to-noise  ratio  at  each  receiver  are  shown  in  the  upper  curve  of 
figure  10.  Receiver  1  utilized  "full"  phase  data  (10  bits  per  phase  sample 
were  actually  used).  Receiver  2  utilized  varying  levels  of  bits  per  phase 
sample  as  indicated  by  the  abscissa  of  figure  10.  The  curve  shows  that, 
for  the  0  dB  input  conditions,  as  phase  quantization  increases  data  transfer 
reduces  faster  than  loss  Increases  until  2  bits  per  sample  is  reached.  Note 
that  for  phase-only  cross-correlation  processing  the  dropping  of  amplitude 
information  at  both  receivers  results  in  a  minimum  loss  of  2  dB  for  the 
0  dB  Input  conditions, 


The  second  data  transfer  reduction  technique  used  half  the  available 
bits  per  complex  sample  for  the  I  component  and  the  other  half  for  the 
Q  component.  Uniform  quantization  bins  symmetrically  spaced  about  0  were 
utilized  with  the  center  of  each  bin  used  as  the  quantizer  output.  To  elimi¬ 
nate  the  use  of  unneeded  quantization  bins  and  to  provide  adjustment  for 
unknown  amplitudes,  the  maximum  magnitudes  of  both  the  I  and  Q  compo¬ 
nents  of  each  receiver's  complex  sample  sequences  were  determined  before 
quantization  and  the  positive  and  negative  quantization  bin  limits  of  each 
receiver  were  adaptively  adjusted  to  the  receiver's  maximum  component 
magnitude.  Results  for  the  case  of  0  dB  input  signal-to-noise  at  esch 
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Figure  11.  Phase  Quantization  Sectors,  2- Bit  Case 


receiver  are  shown  In  figure  10.  Receiver  1  utilized  "full"  I  and  Q  data 
(5  bits  per  I  sample  and  5  bits  per  Q  sample  were  actually  used).  For  the 
0  dB  input  conditions,  data  transfer  reduces  faster  than  loas  increases 
until  2  bits  per  complex  sample  (1  bit  I  and  1  bit  Q)  is  reached.  In  general, 
the  I  and  Q  technique  performs  better  because  it  does  noL  drop  amplitude 
information.  However,  for  the  2-bit  case,  where  minimal  data  transfer  is 
attained  for  a  specified  performance  for  the  0  dB  inputs,  essentially  the 
same  data  are  transferred  from  receiver  2  for  both  techniques.  The  0.  6  dB 
performance  improvement  of  the  I  and  Q  technique  is  attributable  to  the  use 
of  full  I  and  Q  data  from  receiver  1  instead  of  merely  full-phase  data. 

Hybrid  techniques  that  allocate  part  of  the  receiver  1  bits  to  phase  and  part 
to  amplitude  would  provide  intermediate  performance  between  the  I  and  Q 
and  the  phase-only  curves  (e.  g. ,  see  (  8  )  ,  [  9  ]  ,  and  [  10 ) ).  This  I  and  Q 
technique  was  used  because  of  its  simplicity.  It  is  not  optimal.  However, 
as  shown  in  figure  10,  less  than  1  dB  could  be  gained  for  these  conditions 
by  use  of  the  more  complex  data  transfer  reduction  techniques  discussed 
in  the  previous  section.  Note  that  since  the  I  and  C]  technique  provides  less 
performance  loss,  a  specified  level  of  system  performance  would  be  pro¬ 
vided  by  shorter  observation  intervals.  Shorter  intervals,  in  turn,  would 
reduce  the  doppler  decorrelation  loss. 
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CONCLUSIONS 


The  optimal  passive  coherent  location  system  for  the  discrete-time, 
bandpass  case  and  its  performance  characteristics,  including  the  effects 
of  doppler,  were  presented  for  the  fundamental  two-receiver  system. 

These  systems  are  robust  and  make  optimal  use  of  the  observed  information. 
Since  large  quantities  of  data  must  be  transferred  to  implement  these  sys¬ 
tems,  the  performance  losses  incurred  by  two  simple  data  transfer  reduc¬ 
tion  techniques  were  compared  to  performance  loss  bounds  obtained  from 
rate  distortion  theory  and  to  the  losses  of  two  of  the  best  available,  but 
considerably  more  complex,  data  reduction  techniques.  I  and  Q  quantization 
with  simple  uniform,  adaptive  bins  was  shown  to  provide  data  transfer 
reduction  close  to  those  attainable  by  the  more  complex  techniques. 


APPENDIX 


The  Fourier  Transform  of  the  discrete  sequence  of  N  uniformly 
spaced,  complex  samples, 

N 

|  \  =  f(iAt)  | 

i  =  1 
is 

F(u>)  =  ^  fi  e"jwAtl  . 

1  =  1 


Assume  the  f^  samples  are  taken  from  a  bandlimited  process  of  bandwidth 
W  Hertz  and  are  taken  at  or  above  the  Nyquist  rate,  At  <  1/W.  By  inter¬ 
changing  frequency  and  time  in  the  sampling  theorem  for  bandlimited 
signals  [  1]  ,  the  bandlimited  Fourier  transform  may  be  expanded  in  terms 
of  its  frequency  samples  spaced  at  or  above  the  equivalent  Nyquist  rate 
Aw  =  2:r/NAt,  since  the  signal  is  time-limited.  Thus,  the  Fourier  transform 
of  equation  A-l  can  be  expanded  as 
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Fj  =  F(lAw) 


and 


Aw  =  27r/NAt  . 


Consequently,  the  mean  squared  magnitude  of  F(w)  can  be  expressed  as 


N-l  N-l 

,!EIrie  [fi<. 

k=  1  1=1 

k^i 


It  is  well  knowi  that  the  Fourier  frequency  samples,  F^  =  F(i  2?r/NAt),  are 
essentially  uncorrelated  and  their  correlation  decreases  to  0  at  a  rate  of 
at  least  1/N  as  the  number  of  samples,  N,  Increases  ( 11  ] .  Consequently, 
for  moderate  or  large  N,  the  second  summation  of  equation  A-3  will  always 
be  larger  than  the  third  (double)  summation.  Thus,  for  moderate  or 
large  N, 


(A-4) 

The  result  of  equation  A-4  can  be  expressed,  by  the  use  of  equation  A-l,  as 
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